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Abstract 

A quasi-static process is realized in a purely quantum-mechanical model which 
is described by oscillator (or particle) systems having relative-phase interactions. 
Time development of a mixture of two oscillator (or particle) systems which have 
different temperatures is investigated in the quantum mechanical framework. We 
study how they go to a thermal equilibrium. Time dependence of the temperature of 
the object is numerically studied and the difference of the time dependence between 
classical objects and quantum objects is clearly pointed out. 
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1 Introduction 



In the papers published recently [1,2] we have proposed a quantum mechanical description 
of thermal equilibrium states in terms of eigenstates of relative-phase interactions H}([9]) 
between oscillators (or particles) with phases which are introduced as aj{9) — aje~'^^\ 
where aj is the annihilation operator of the jth oscillator (or particle). The interaction 
projects the physical space TCq of the original Hamiltonian Hq = ^^f^iOjCij, which is 
spanned by the direct products of the usual number states Y[f=i 1% onto the thermal 
subspace Tithermai of TYq, whcrc all the eigenstates having the same energy eigenvalue have 
the same probability as if the fundamental principle of thermal equilibrium, the principle 
of equal a priori probability, is satisfied. It has been shown that the thermal subspaces 
characterized by the different phases are equivalent as the thermal subspace of TIq [3] . The 
microcanonical ensemble is derived by the average operation with respect to the whole of 
the equivalent subspaces which cover the original physical space TYq completely. In the 
model it has been shown that the effective density-matrix for arbitrary small subsystems 
of the total system are approximated by canonical ensemble of usual statistical mechanics 
[4]. In the procedure a common temperature is uniquely introduced and determined by 
the physical quantities of the heat-bath which is taken as the remaining part of the total 
system except the small subsystems. We can also see that this model can be described by 
using new operators defined on the thermal subspace V-thermai in very simple forms. Fur- 
thermore in this description we can easily see that the relative-phase interaction provides 
transitions only on the thermal subspace. This means that quasi-static processes [4] can 
be defined as quantum processes. In this paper we shall show the above results and apply 
it to the calculation of change of temperatures. 

In §2 we briefiy review the model and derive useful operators on the thermal subspace. 
The interactions of two systems and quasi-static processes are investigated in §3. In §4 
time evolutions of the two systems are discussed. Time evolutions of the temperatures of 
the interacting systems are calculated in §5. Remarks on the model are presented in §6. 



2 Basic ideas of the model and useful operators on 

'^thermal 

In order to make our discussions clarify we review on the previous work describing thermal 
equilibrium for an oscillator system with one energy scale e. (For details, see Ref. 1.) We 
shall also introduce useful operators on the thermal subspace. 

2-1. Model 

The model is described by the system composed of oscillator (or N particle) which 
are described by the Hamiltonian Ti^ = e X^^i ojcj, where and Oj arc, respectively, the 

creation and annihilation operators and follow the commutation relations [aj,a]] = 
The phases are introduced in terms of the replacement of the operators with the new 
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ones aj{9) — aje and 0^(9) — aijC^^^ . In such a replacement Hq and the commutation 
relations do not change. Using the usual number states satisfying OjUjlnj >= nj\nj > 
and aj\0 >= 0, we can obtain the eigenstates of Hq with the eigenvalues Em — 
(M = 0,1,2,...) as follows; 

TV 

|M; [nj], [^j] >= |M;ni,n2,--- , njv; ^i, ^2, • • • , ^Jv >= H '^J' ^ ^ELi^fc.M > (1) 

where \nj,6j >= \nj > e"*^^^ . Note that wc have the equation aj(^^)|nj, >= y^nj~+T|nj+ 
l,9j >, aj{9)\nj, 9j >= ^Jfrj\nj — l,9j > and aj{9)aj{9)\nj,9j >= nj\nj,9 >. States in 
which all the eigenstates with the same energy eigenvalue have the same probability, as if 
the principle of equal a priori probability with respect to the eigenstates of Hq is satisfied, 
are given as 

|M, TV, [9] >= J2 \M; [9,] > /^W{M,N), (2) 

p[nj] 

where the sum should be taken over all the different combinations of [nj] ~ {ni,n2, tt-at) 
and the number of the combinations W{M, N) is given by W{M, N) = ^^H^N-iy ■ 

2-2. Relative-phase interaction 

We introduce the following operators [1-3] ; 



a]{9) ^ al{9){^ N, + a,{9) ^ + l)-\{9) for j = 1, TV, (3) 

where Nj = a}-{9)aj{9) is the excitation- number operator for the jth oscillator. We see 
that 

a]{9)\nj, 9j > = {rij + 1, 9j >, for Vn^- = 0, 1, 2, ... 

aj{9)\nj, 9j > = {rij - 1, 9j >, for Vn^- > 1 (4) 
and aj{9)\0 >j— 0. The relative- phase interaction with the energy scale e^is written as 

/V ^ 1 ^ 

where a] = Q;](0), = 0:^(0), iV = J2f=i-^j (^^^ total excitation-number operator) and 
the relative phase 9jk — 9j — 9k- We obtain the relation 

H\M, N, [9] >= (e + eg)M\M, N, [9] > (6) 

for H = Hq -\- Hj{[9]). We see that the physical space of the total Hamiltonian H is 
spanned by the set of the states \M,N,[9] > for M = 0, 1, 2, which is nothing but the 
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thermal subspace Hthermai of the original space Hq. Thermal limit will be described by 
the limit €g 0, where the eigenvalues of H can approximate to those of Hq. Note here 
that the states with the different phase parameters \M,N, [6'] > are not the eigenstates 
of H}{[9]), that is, H}{[e])\M,N, [9'] egM\M,N, [9'] > for [9'] 7^ [9]. 

2-3. Derivation of canonical ensemble and temperature 

We divide the total system into two groups, hj = aj (j = 1,2, ...,Nh)) (6-group) and 
dj^N^ = ttj for j > Nh + l ((i-group). The eigenstates of the total Hamiltonian |M, N, [9] > 
can be written in terms of the direct product of the eigenstates of the groups such that 

M M / w( M — AM 

® ^W{M,,N,)\Ma, Nd, [0] >d, (7) 

where \Mb, N^, [9] >{, and \Md, N^, [(/)] >d are, respectively, the eigenstates on the thermal 
subspace of the 6-group and that of the rf-group. Hereafter we shall treat the 6-group as 
the heat bath and the (i-group as the object. Then the 6-group is taken to be much lager 
than the d-group, i.e., Ni, » Na and Mj, » M^. 

In order to evaluate the density matrix for the object (d-group) we have to perform the 
partial trace operation (Tr/), because any quantum numbers of the heat-bath (6-group) 
state are not measured in this processes. We obtain the effective density matrix for the 
object as follows; 



^ ^M,+Ma,M m \Md:Nd, [(p\ >d d < Md,Nd, [(p\\. 



PeffM) = Trj{\M,N, [9M X M,N, [^][0]|), 

WiM, N) 

(8) 

In the limit where the interaction with the heat bath is small enough to ignore, i.e., 
ep << e, we may again introduce the average over the un-measurable phases of the object, 
[0] , because H^^ — e '}2k=i ^l^^^ [0]-dependence. Thus we obtain the density matrix 

for the mixed states with respect to the eigenstates of H^^ 

^ ^ W( M N ) 

Peff = E E ^M,+Mrf,M E l^'^; K] >d d < Md] K]|. (9) 

Note here that in the above procedure the average operation with respect to the phases of 
the heat bath, [9], are not carried out, while it must be done to derive microcanonocal en- 
semble for the heat bath as was discussed in Refs. 1-3. The canonical ensemble is realized 
at the maximum of the probability Pm{Mi,) = W{Mi,, Nj,)W{M -Mb,N- N^) /W{M, N) 
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for finding the state jM^, N^, [9] >b <8)|Md, A^^, [0] >d in |M, N, [6'][0] > [4]. At the maxi- 
mum we can estimate the change of the coefficients in the probabihty with respect to the 
change of the object energy, that is, AE"'' = eAM^ with AMa Ma « M, Mb, as 

W{Mb- AMd,Nb) ^ Mb ^ -lAE'X'Hi+N./M,) /-i^N 

W{Mb,Nb) ^M^b + Nb' ■ ^ ^ 

Thus we can determine the temperature T as 

kT^€{\n{l + Nb/Mb))-\ (11) 

It is obvious that the temperature is written down by using the quantities of the heat 
bath, that is, the mean excitation number (Mb/Nb) of the heat bath. We obtain the 
density matrix for canonical ensemble in the limit M ~ Mb, N Nb ^ oo 



oo 

Tob 



Md=0 p[mj] 

As for the details of the discussions presented in this section, please see §5 of Ref. 1. 

2-4.- Useful operators on the thermal subspace 

Here we shall present a pair of useful operators in the following discussions, which are 
defined by 

N N 



Ae = \In + nY,o^^{0)IN, Ae = a,(^)(\/i^+7V)-\ (13) 



where aj{6) = aj{9)y Nj + 1. It should be noted that they satisfy the commutation 
relation [^e,^^] = 1. We easily se that on the thermal subspace Hthermai 



Ae\M, N, [9] > = VM\M - 1, A^, [^] > 



Ae\M, N, [9] > = VmTT\M + 1, N, [9] > (14) 

and AeAe\M, N, [9] >— M\M, N, [9] >. We may say that they each other behave as her- 
mitian conjugate operators on 'Hthermai- Now we can write the relative- phase interaction 
and the eigenstates in very simple forms on Hthermai as follows; 

= ^gAeAe, \M, N, [9] >= -^{Ae)''\0 > . (15) 

An important feature of the new operators can be seen in the choice oi N — 1, where 
they coincide with the original operators such that Ae — a{9) and Ae — a^{9). 
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3 Interaction between two systems and quasi-static 
process 

Let us consider interactions between two oscillator systems composed of the same os- 
cillators (or particles). In order to distinguish them one is called a-system described 
by the original Hamiltonian Hq = e X^^i OjCij and the other 6-system described by 

Hq — ^Yl!j=i^\^r Hereafter we shall represent the operators defined in 2-4 by Ae for 
the a-system and by for the 6-system. We postulate that the a-system is in the ther- 
mal equilibrium on 'Hthemmi 6-system in that on Ti-^thermai before the two systems 
are mixed. That is to say, the a- and 6-systems are, respectively, described by the states 
\Ma,NaAO] >= T^W'^IO >a and \M,,N,,[<P] >= -^{B^)^>'\Q The thermal 
interaction may be written as follows; 

Hi = eg At At, (16) 

Where At = Ae+B^, At = Ae + B^. The terms Hf{[9]) = e^AeAe and H]{\(\)\) = egB^B^ 
involved in Hj are, respectively, nothing but the relative-phase interactions for the a- 
and b- systems, which are already taken into account in the description of the states 
\Ma, Na, [9] > and \Mb,Nb,[(j)] >. Then the interaction between the a- abd b- systems is 
represented by 

Hf = eg{AeB^ + B^Ae). (17) 

From the relations of (14) we see that this interaction induces the transitions only among 
the states on the direct-product space of two thermal subspaces for the two systems such 
as |M„, Na, [9] > ®|Mb, Nf,, [0] > on Hthermai ® ^Lrma^ ^ud then it preserves both of 
the thermal equilibrium for the a- and fo-systems, even though the thermal equilibrium of 
the total system is not realized. It means that the processes induced by this interaction 
represent quasi-static processes. Then we may say that all thermodynamical quantities of 
the two systems, such as the temperatures of a- and fo-systems, chemical potentials and so 
on, are well-defined physical quantities at arbitrary moments of this physically changing 
process. 

4 Time evolutions of the two systems 

In the interaction picture the unitary time evolution operator is evaluated by the formula 

U{t,0) = exp[-i [ Hint{t')dt'/h] (18) 
Jo 

where 

= €g{Ae{t)B^{t) + B^{t)Ae{t)), (19) 
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where = + Hf{[d]), Aeit) = Aee-'^'+^s^*/^,A0it) = ^^e^(^+^«)*/'^ and so on. 
Note here that the t-dependence drops from Hint{t) because the operators Ag{t) and 
B^{t) have the same time dependence and then U{t,0) is the exact formula because 
[Hint{t'), Hintit")] = is satisfied for arbitrary t' and t" . Putting ujg = Cg/h, we obtain 

U(t, 0) = eM-iugiAeB^, + B^A)^]- (20) 

We study the case where two systems are mixed at t = and the initial state is given by 

>= \Ma, Na, [6] >0\Mb, Nb, [0] > . 

We should study two different situations; one is expressed by the interaction between 
two systems which have a comparable size, that is, A^a ~ A^^, and the other is the case 
where one of the two (heat bath) is extraordinary lager than the other (object), e.g., 

Na » Nh, Ma » Mh. 

Case (1) for Na ~ Nh'. We may consider that all the constituents in the both systems 
couple each other. The density matrix is written down as p = U{t, 0)|\1'/ > < ^/|t/(t, 0)^" 
and the effective density matrix for the object (6-system) is obtained by performing the 
partial trace operation (TV/) for the variables of the other (a-system) as 

M 

Peff = Yl K{t)\''\m,Nb, [4>] > < m,Nb, [0]|, (21) 

m=0 

where the coefficients Cm{t) are defined by U{t,0)\^i >= ^^=o'^m(^)l^ ~ itt-, Na, [9] > 
®\m,Nb, [0] > with the constraints I]m=o |cmWI^ = 1 and M ^ Ma + Mb. The effective 
density matrix pe// is represented by a mixed state on 'H^hermai- We should pay attention to 
the fact that the partial trace operation with respect to the variables of the non-objective 
matters (a-system) plays an essential role in deriving the mixed state property of peff- It 
is also note that Pe// is not yet represented by any mixed states on T-Lq. The decoherence 
with respect to the eigenstates of H^, that is, the decoherence on TYg, is accomphshed by 
the average over the phases [0] as noted in 2-3. In the following discussions, however, 
we have to use pe// because we cannot take the thermal limit eg 0, where any time 
evolution does not take place. 

We can study the time development of the total system in perturbation because of 
e» Eg. The matrix elements &X, t — t are obtained as follows; 

< -^ipir, 0)1*, > ~ 1 - ]^{ugTf{2MaMb + M„ + M^), 
<Ma + l\ <Mb- l|C/(r,0)|*/ > ~ -iujgT^iMa + 1)M„, 

<Ma-l\ <Mb+ 1\U{t,0)\^i > ~ -iujgT^{Mb + l)Ma (22) 
and others are less than the order of {(jgr)"^, where \Ma > \Mi, >= \Ma,Na,[9] > 
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(8)|M6, Nb, [(f)] > is used. The probabilities up to the order 0{{u;gTy) are given by 

P{Ma, Mb) ^ 1 - iuJT)\2MaMb + M„ + Mb), 

P{Ma + l,Mb-l)- Kr)2(A4 + l)Mb, 

P{Ma - 1, M5 + 1) ~ Kr)'(M6 + 1)M„. (23) 

We obtain the time dependence of the mean excitation number of the object as 

< B^B^ >= Mb(l) ~ Mfc - Kt)2(M6 - M„). (24) 

This result is quite reasonable because the energy of the object increases if it is smaller 
than that of the environment and vice versa. Furthermore the change of the energy 
stops when both of the energies become same. Let us estimate the time when the time 
evolution stops. The time evolution equation for the excitation numbers is obtained from 
the relations for the nth steps of the perturbation 

Mfe(n) = Mb{n - 1) - {LUgT)^[Mb{n - 1) - M,(n - 1)], 

Main) = Main - 1) - iu;gT)^[Main - 1) - M^in - 1)]. (25) 

We have the equations for the difference m(n) = M^in) — Main) as follows; 

Am(n - 1) = -2iugT)'^min - l)An, (26) 

where Am(n — 1) = m(n) — m(n — 1) and An = 1. Since we may treat n as a continuous 
variable in the region n >> 1, we derive the relation 

m(n) ~ (Mb - M„)e-^'"^*", (27) 

where the time corresponding to the nth step t„ = rn. Prom the above result and the 
relation for the energy conservation Main) + M^in) — Ma + Mb for Vn, we get 

Ma + Mb Mb -Ma o,,2^, , , 

Mbitn) = \ + ^ ° e-^"^"^". (28) 

The time when the thermal equilibrium is almost realized is that when the difference 
m(n)/Mb approximates to 0. Then we can estimate it from the relation m(n) ~ 1 << 
Ma, Mb as 

titmal^^H\Mb-Ma\). (29) 

Note that the exponential damping property of Mbitn) for tn is same as that of classical 
damped oscillation. We may put ujgT < (M^ + M^)"^/^ because the consistency of the 
perturbation, that is, all the probabilities in (23) must be positive. We have 
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Prom (11) the temperature of the object after the n-perturbation steps is derived as 

T,(n) ~ (e/fc)[ln(l + 2N,e/{E, + E, - {E, - Eje-^-^*"))]"!. (31) 

Since the time {tmeas) required for measuring the energy with the accuracy Ai? ~ e is 
much shorter than the perturbation time (tpert), i-c, tmeas ~ ^/e << tpert ~ ^/^g because 
of e >> tg, we have time enough to measure the temperature which is determined by the 
knowledge on Hli^^^^ai involved in Hq. 

Case (2) where an object is put in a huge heat hath: We have to take account of the 
fact that the object interacts with a very small part of the heat bath. Therefore the heat 
bath must be divided into two parts and the part (a-system) which interacts with the 
object (6-system) is comparable to the size of the object, that is, Na — N^, and very much 
smaller than the non- interacting part. Such a separation has already been carried out in 
2-3 and written in the canonical ensemble given in (12). The effective density matrix for 
the object can be obtained as follows; 

OO Ma+Mt 

Peff = E E e-f'(''^+^^-^>\Ut)\'\m, N,, [0] >< m, AT, < [0]|, (32) 

Ma=0 m=0 

where ,5 = 1/ kT (T is the temperature of the heat bath). Following the same pcrturbative 
treatment, we derive the change of the mean excitation number of the object as 

< B^B^ >~ M, - {ugTf{M^- < M„ >), (33) 

where the expectation value for the excitation of the interacting part is evaluated as 

<Ma>-Y. - (34) 

The difference of the result of (33) from that of (24) appears in the fact that < > in 
(33) does not change during the interaction with the object, while in (24) does. The 
change of the excitation number after the n-perturbative steps is obtained as follows; 

Mt(n) =< Ma > +{M,- < Ma >)((1 - (1 - (a;,r)2)«). (35) 

In the limit of n — > oo we see that 

\imt^^^Mb(n) ^<Ma> . (36) 

It means that the object goes to the thermal equihbrium provided by the heat bath. Since 
we may represent Mh{n) by Mi,{tn) ^< Ma > +{Mh- < Ma >)e~'^9'^*" for large n, the 
time when the thermal equilibrium is almost realized is estimated as 

tfLmal ^ ;lzH\M,- <Ma> |), (37) 
9 
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which is just twice of the (!)■ This is quite reasonable because in this case 

the a-system is a part of the heat bath does not change the thermal properties, whereas 
the a-system changes the states towards the thermal equilibrium in case (1). That is to 
say, the change of the temperature of the 6-system in this case must be twice as large as 
that of case (1). We shall easily be able to observe this difference in experiments. 

Note here that Mb(tn) damps exponentially with the time for t„ >> l/^^gT hke the 
classical damped oscillation, while it damps as Gaussian for tn << l/ci^g. The perturbation 
approach used in this section is one of the weak coupling limits [5] . 



5 Some numerical considerations 

Here let us numerically study the time-dependence of the temperature which may be the 
most probable observation in experiments. We investigate a simple example for case (1), 
that is, we study the change of the temperature for the mixture of two same particle 
systems (gases) which have a similar weight {Na — Ni,) but different temperatures {Ma 7^ 
Mij) and bottled in the same-sized cubes (the volume=/^ (cm^)). Since the energy spectrum 

of a constituent particle are derived as En^^ny,n^ = + + ^2); where m is the mass 

of the particle and Ux, riy, Uz are positive integers, all the energy levels are represented 
by the lowest energy [Eq) x integers and then we may take Eq as the fundamental level 
spacing e used in the previous discussions. (Exactly speaking, we have to use the method 
for the model involving different energy scales presented in Ref. 2.) The lowest energy is 
obtained as 

e = £;o~0.25x 10-'°^(erg), (38) 
where the gas is taken as the hydrogen gas. The temperature is written down as 

where Mft(t) is given by (28) and t — mis used. 

Let us study a case corresponding to our daily phenomena, say I — 1cm and T ~ 
10^ {°K). Since the factor 0.2 x 10~^^ is very small, the approximation \n{l-\-Ni,/Mi,(t)) ~ 
Nh/Mh{t) « 1 must be realized even at very low temperatures. Thus we get 

T(t) ^l{Ta + n) + 1(T, - T,)e-2(-^)* (40) 

for t » l/UgT, where and T5 are, respectively, the initial temperatures of the a- and 
6-systems. This formula is quite reasonable from the classical point of view. We can 
determine the parameter cUgT from experimental data. 

Let us study the case for the quantum size, for instance, I — lO^^cm. We get 
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The difference from the above result (40) is found only below the critical temperature 
7; ~ 20/ln2 (~ 29) where the relation N^/M^{t) » 1 must be satisfied. We have 

nit) ~ Tdf) + — X rJ/)^ ^'~^° e-^("^)^ (42) 

for t » 1/ujgT, where Th{f) = e/[k hi{{Na + Nij)/{Ma + Mb))] representing the tempera- 
ture in the final thermal equilibrium may approximate to 20/ln(2iVb/(Ma + M^)) because 
of Na ~ N),. We see that the time-dependence is same as that of the classical case but 
the factor of the time-dependent term is quite different from that of the classical case. 
Note that the numeral 1/20 in (42) has the size (/(cm))dependence such as 5 x 10^"^/^. 
The difference of the quantum case from the classical case can experimentally be seen 
by measuring the size dependence. For case (2) we can proceed the same argument by 
replacing Ma with < Ma >. 



6 Remarks 

We have showed that a simple non-thermal equilibrium phenomena can be treated in 
the present quantum mechanical framework for thermodynamics. It should be stressed 
that in this scheme we can proceed all evaluations in preserving quasi-static property. 
We can, therefore, determine temperatures of all subsystems at arbitrary time. Davies 
presented quantum theory of open systems [5]. In his approach, however, we have no 
definite way to determine the temperature of the object at arbitrary time. It is also note 
that in his approach the canonical distribution for the heat bath is put by hand, whereas 
the distribution is derived in the large heat-bath limit (M and 00) in the present 

scheme as was shown in (12). Exactly speaking, we can study the deviation from the 
canonical distribution for systems having not large M and A^. It is interesting that our 
average operation over the phases has the same role as that of the projection operator Pq 
introduced in section 10.2 of Ref. 5, which projects the space of the closed system onto the 
subspace described by real diagonal matrices. Thus we see that the phases in our model 
play an essential role for the derivation of the microcanonical and canonical ensembles. 
(See §2-3.) There has been no model in which temperatures of interacing systems are 
quantum mechanically calculable. In fact we need not to introduce any temperatures a 
priori in our scheme. This point is quite different from other approach, for instance, the 
model for an oscillator interacting with environment by Unruh and Zurek [6]. 

The magnitude of the relative-phase interaction is still unknown. One possible obser- 
vation of the interaction in cold- neutron experiments has been proposed [3] . We can also 
estimate the strength of the relative-phase interaction by measuring UgT through the time 
development of temperatures. Whether its magnitude is universal or not will be a good 
test for the universality of the relative-phase interaction. We should, however, not forget 
about the fact that in the present scheme the viscosity introduced from interactions on 
surfaces (boundary of heat bath and object), which will disturb the thermal equilibrium 
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of the object, is not taken into account. Those problems must be studied in more general 
formalism [5]. 
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